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In this paper we report for the first time on the necessity of the refinement of the concept of 
generahzed chaotic synchronization. We show that the state vectors of the interacting chaotic 
systems being in the generahzed synchronization regime are related with each other by the functional, 
but not the functional relation as it was assumed until now. We propose the phase tube approach 
explaining the essence of generalized synchronization and allowing the detection and the study of 
this regime in many relevant physical circumstances. The finding discussed in this Report gives a 
strong potential for new applications. 
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Chaotic synchronization is one of the fundamental phe- 
nomena, widely studied recently, having both theoretical 
and applied significance [1]. One of the interesting and 
intricate types of the synchronous behavior of unidirec- 
tionally coupled chaotic oscillators is generalized synchro- 
nization (GS) [H,!!!. This kind of synchronous behavior is 
said to mean the presence of a functional relation between 
the drive and response oscillator states [3, HI and has been 
observed in many systems both numerically [^-[11 and ex- 
perimentally [ol-llH, with many interesting features (7l,fT2| 
and possible applications [l3|, 141 of this regime being re- 
vealed. 

The definition of the GS regime generally accepted 
hitherto is the presence of a functional relation 



yW = F[x(i)] 



(1) 



between the drive x(t) and response y{t) oscillator 
states 0, Having based on this definition the dif- 
ferent techniques for detecting the presence of GS be- 
tween chaotic oscillators had been proposed, such as the 
nearest neighbor method 0, EH , the auxiliary system ap- 
proach H or the conditional Lyapunov exponent calcu- 
lation [5| , with the auxiliary system approach being gen- 
erally the most easy, clear and powerful tool to study the 
GS regime in the model systems, whereas for the analy- 
sis of the observed experimental time series the nearest 
neighbor method, as a rule, is more applicable [11]. 

In this Report we report for the first time on the ne- 
cessity of reconsidering and refining the existing concept 
of generalized chaotic synchronization. The main rea- 
son of this refinement is the following. Let x(to) = xq 
and y(to) = yo be the reference points belonging to the 
chaotic attractors of the drive and response oscillators 
being in the GS regime, respectively. For the neigh- 
bor point x(ti) = Xi of the drive oscillator such that 
I |x^ — xo| I < £ its image y{ti) = in the response sys- 
tem is also close to the reference point yo (see Q for 
detail), i.e., ||yi — yo|| < ^{s)- Having linearized Eq. ([T]), 
one obtains that 



where J is the Jacobian operator. Since the form of 
the functional relation F[-] can not be found explicitly 
in most cases, Eq. ([2]) may be rewritten in the form 



(3) 



where A = JF[xo] is the unknown matrix and 
^x^ = Xi — xo , 6yi = Yi — yo are the vectors characteriz- 
ing the deviation of the points under consideration x^, 
from the reference points xq and yo, respectively. With- 
out the lack of generality we shall suppose below the 
identical dimension m of the phase space of the drive 
and response systems. 

Although the coefficients of the matrix A are unknown, 
the validity of Eq. ([3|) may be verified if there are N > m 
nearest neighbors x^ of the reference point xq and corre- 
sponding them vectors y^ of the response system. Hav- 
ing tested the presence of the generalized synchronization 
(e.g., with the help of the auxiliary system approach) we 
can pick out m nearest neighbors x^ {i = 1, . . . , m) and 
corresponding to them vectors y^ to determine the co- 
efficients aij of the matrix A with the help of Eq. ([3]). 
To reduce the influence of the inaccuracy we have to se- 
lect such vectors x^ (and ^x^ = {Sxn, . . . , 6xim)^ ^ respec- 
tively) from the whole set of N vectors for which 



I det(X)| = max. 



where 



X = 



SX21 



Sxi2 
SX22 



, . . Sxim \ 
■ . . Sx2m 



(4) 



(5) 



y SXfYil ^XfYi2 ' • • ^XfYim J 



Having determined the matrix A we can now find the 
vectors Szi, (i = m + 1, . . . , TV) as 



6zi = A(5xi, 



(6) 



Yi - yo = JF[xo](xi - xo), 



(2) 



and compare them with the vectors Sji of the response 
system (or compare vectors = yo + Szi with y^) to val- 
idate the correctness of Eq. (j3|). 
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FIG. 1: (Color online) The vectors yi (■) and (o) of the 

response Rossler system dT)) for £ = 0.3. The numbers i of the 
vectors are shown by the regular and italic fonts, respectively 



Altought, at first sight, it seems that there are no fun- 
damental causes due to which Eq. (j3j) may fail, in reality 
Eq. (|3]) is not correct. To illustrate this fact we have stud- 
ied numerically the synchronous behavior of two coupled 
chaotic Rossler oscillators 

Xd = -^dVd - Zd, Xr = -UJrVr - + ^{Xd - Xr), 
yd = UJdXd + Ciyd, Vr = ^rXr + a^r, 

Zd=P^ Zd{Xd - C), Zr = Zr{Xr - c), 

(7) 

where x = {xd^Vdi ^d)^ [y = (^r, ^r, ^r)"^] are the carte- 
sian coordinates of the drive [the response] oscillator, dots 
stand for temporal derivatives, and £ is a parameter rul- 
ing the coupling strength. The other control parameters 
of Eq. ([7]) have been set to a = 0.15, p = 0.2, c = 10.0, 
in analogy with our previous studies [sl, [l6| . The cj^- 
parameter (representing the natural frequency of the re- 
sponse system) has been selected to be ujr = 0.95; the 
analogous parameter for the drive system has been fixed 
to ujd = 0.99. For such a choice of the parameter values 
the boundary of the generalized synchronization regime 
found with the help of the auxiliary system approach is 
around sgs ^ 0.11. 

Having chosen the reference point xq of chaotic attrac- 
tor of the drive oscillator randomly, one can find its near- 
est neighbors x^ {i = 1,...,A^) (and corresponding to 
them vectors of the response system), select (accord- 
ing to Eqs. dH and (|5|)) the vector basis xi_3 to determine 
the matrix A and check condition (|3]) with the help of 
Eq. ([6]) and the rest of vectors x^, y^, (z = 4, . . . , N). 

In Fig. [1] the vectors (i = 4, . . . , 10) obtained with 
the help of Eq. (|6|) as well as the vectors y^ of the response 
system are shown for the coupling strength s = 0.3. 
The value of the coupling strength exceeds greatly the 
threshold ecs of the generalized synchronization, the 
GS regime demonstrate great stability, and, as a conse- 
quence, Eq. (j3j) is expected to be correct. However, con- 
trary to expectations, the vectors and y^ differ from 
each other sufficiently testifying that Eq. ([3]) fails. As 
a matter of fact, the failure of Eq. (|3|) is also observed 




FIG. 2: The schematic representation of the nearest vectors 
Xi(t), the phase trajectories Xi(s) and the phase tube Tr(t) 

for other reference points of the drive Rossler oscilla- 
tor as well as for other chaotic dynamical systems (e.g., 
Lorenz oscillators). Since Eq. (|3]) is just the linearization 
of Eq. ([1]), the failure of Eq. ([3]) is the evidence of the 
incorrectness of Eq. ([1]) being the main definition of the 
generalized synchronization concept. At the same time, 
plenty of results obtained hitherto are in the very good 
agreement with the generally accepted concept of GS. It 
means that the concept proposed by N. Rulkov et al. Q 
works in some circumstances, but, in general, must be 
refined. 

The core idea of this correction is the following. The 
state of the response system y{t) depends not only on the 
state of the drive oscillator x(t) at the moment of time 

but on the history of the evolution of the drive system 
during time interval (t — r, t] as well. Indeed, according 
to the concept of GS, synchronization means that the re- 
sponse oscillator y{t) comes to the state defined uniquely 
by the drive system, with the convergence time r being 
connected with the largest conditional Lyapunov expo- 
nent A^, i.e. T ^ l/\Xi\. In other words, F[-] in Eq. ^ 
must be considered as a functional^ but not a functional 
relation. Obviously, in this case Eq. ([3|) obtained un- 
der assumptions that F[-] is the functional relation is not 
satisfied as it has been shown above (see Fig. [1]). 

Considering F[-] as the functional^ one have to replace 
Eq. ([2D by 

t 

Sy^{t)= J JF[Ms)]SMs)ds. (8) 

t-T 

Having supposed that the deviation (5x^(s) from the ref- 
erence trajectory xo(s) {t — r < s < t) is small, in view 
of the linearity one can write 

(5X^(5) = B(5)JXi(t), t-T<S<t, (9) 

(where B(s) is the matrix with the time-dependent coef- 
ficients) that results in 

t 

Syi{t)= J JF[Ms)]B{s)S^i{t)ds. (10) 

t-T 
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ll^y^^) ~ ^2*(^))ll/ll^y*(^)ll (^) ^he asynchronous dynam- 

FIG. 3: (Color online) The vectors y. (■) and (o) of the ^ = ^'^^ ^""^ generalized synchronization regime, 

1 X /TtI r n o 1 f £ = 0.3. The histoCTams havG been obtained for the resDonsG 

response Kossier system (|7|) tor e = 0.3, the length or the ^ —.^ ^ ^ i r i i i • ^r^r^ 

1 . 1 • ir>r> rpu u • r ^ Kossler system (ITD, the length or the phase tube rs r = 100 

phase tube rs r = 100. Ihe numbers i ot the vectors are j t> 

shown by the regular and italic fonts, respectively 



and, as a consequence, in 

6y,{t) = C{t)S^i{t), (11) 

where C{t) is the square (m x m)-matrix defined as 

t 

C{t)= J JF[^i{s)]B{s)ds. (12) 

t-r 

So, Eq. (fTTj) coincides formally with Eq. (j3j) and, there- 
fore, it may be also validated by the calculations of 
vectors in the same way as it has been done for 
Eq. (|3]). At the same time, Eq. (|3]) has been obtained 
under assumptions that the vectors xo(t) and x^(t) are 
close to each other, whereas Eq. (pT|) has been obtained 
for more stricter restriction requiring the nearness of 
the trajectories xo(5) and Xi{s) during the time interval 
t — r < s < t. Since for the chaotic systems the phase 
trajectories can converge in one direction of the phase 
space and diverge in another one, the neighbor vectors 
xo(t) and x^(t) may be characterized by the very distinct 
phase trajectories xo{s) and Xi{s) for t — r < s < t. The 
schematic representation of such a situation is given in 
Fig. [21 Although the vectors xi(t) and X2(t) are close 
to the reference point xo(t), only the vector X2(t) obeys 
Eq. (pTj) due to the nearness of the phase trajectories 
xo(5) and X2(s), whereas for the vector xi(t) Eq. (pTj) 
fails, since the phase trajectory xi(s) is not close to 
the reference one xo(s) during the whole time interval 
t — r < s < t. Therefore, to verify Eq. (pT]) we have to 
consider not all vectors Xi{t) being nearest to the refer- 
ence point xo(t), but only vectors which are characterized 
by the phase trajectories x^(s) being close to the refer- 
ence one xnfg). Having based on the idea of phase space 
strands [13, [l8| , to eliminate the ineligible vectors (like 
xi (t) in Fig. [2]) we introduce into consideration the phase 
tube 

Tr{t) = {x : \xoj{s) - Xj \ < 6^,17^1 ,5 G [t - r; t]} (13) 

and take into account only vectors whose phase trajecto- 
ries pass through this phase tube (like X2(t) in Fig. [2]). 



The result of this examination for Rossler systems ([7]) 
with the same set of the control parameter values and 
the coupling strength as before is given in Fig. [3l the 
length of the phase tube is r = 100. One can see that 
the calculated vectors (t) are in the excellent agreement 
with the vectors (t) of the response Rossler system that 
confirms both the correctness of Eq. (pTj) and, as a con- 
sequence, the statement that F[-] is the functional^ but 
not the functional relation. 

With the increase of the coupling strength between 
chaotic oscillators the absolute value of the largest con- 
ditional Lyapunov exponent grows, whereas the time 
interval r (the length of the phase tube Tr{t)) decreases. 
Finally, in the lag synchronization (LS) and complete 
synchronization (CS) regimes the value of r tends to be 
zero. Therefore, in the LS and CS regimes Eq. ([3]) is 
satisfied for all neighbor vectors x^(t) without any addi- 
tional requirements concerning the phase trajectory near- 
ness. In other words, the state vectors of any chaotic 
systems being in the GS regime (but not in the LS or 
CS regime) are connected with each other with the func- 
tional, whereas in the LS and CS regimes (which are the 
strong form of GS) they are related with each other by 
the functional relation. 

Though the phase tube approach has been here applied 
to the model systems, we expect that it can be used in 
many other relevant circumstances. Since the statistics 
for the difference between Syi{t) and Szi{t) vectors are 
radically different for the synchronous and asynchronous 
motion (see Fig. 2]), the important feature of this ap- 
proach is the possibility to consider the relation between 
vectors (pTj) for the analysis of the registered experimen- 
tal data (vector or scalar, using the Takens approach [l9|) 
when the other classical methods of GS detection are in- 
accurate or unapplicable. Moreover, the proposed ap- 
proach may be used as the method to detect the GS 
regime, including the case when the chaotic oscillators 
are coupled mutually, since all arguments given above are 
also applicable for the case of the bidirectional coupling. 

To prove the generality of our findings we have also 
studied numerically two mutually coupled generators 
with tunnel diodes ^j. In the dimensionless form the 
dynamics of such generators is described by the equa- 
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FIG. 5: (Color online) The vectors yi (■) and (o) of the 

second generator with tunnel diode ([14]) for e = 0.15. The 
numbers i of the vectors are shown by the regular and italic 
fonts, respectively, (a) All neighbor vectors are used; (6) only 
vectors whose phase trajectories pass through the phase tube 
with length r = 110 are used 

tions dO, [in 

^1,2 = ^i,2[^(^i,2 - £{y2,i - Via)) + Via - ^1,2], 

yi,2 = -^1,2 + £{y2,i - ^1,2), (14) 

Mil, 2 = ^1,2 - /(^i,2), 

where /(O = + 0.002sh(5^ - 7.5) + 2.9 is the dimen- 
sionless characteristics of nonlinear converter, h = 0.2, 
/i = 0.1, uJi = 1.09, Lj02 = 1.02 are the control parameter 
values, e is the coupling parameter strength. The indexes 
"1" and "2" correspond to the first and second coupled 
systems, respectively. For such values of the control pa- 
rameters the threshold of the generalized synchroniza- 
tion regime determined by the moment of the transition 
of the second positive Lyapunov exponent in the field of 
the negative values [lH, [23| is around sgs ^ 0.08. 

As in the case of Rossler systems considered above we 
have chosen the reference point xq of chaotic attractor of 
the first oscillator randomly and analyze the behavior of 
its nearest neighbors (i = 1, . . . , A/") and corresponding 
to them vectors and . The choice of the vector basis 
xi_3 has been performed in the same way as in the case 
considered above. 



In Fig. [5] the vectors (■) and (o) of the second 
generator with tunnel diode (p!4|) for the coupling pa- 
rameter strength s = 0.15 exceeding greatly the thresh- 
old value of the generalized synchronization regime onset 
egs are shown. Fig. corresponds to the case when 
all neighbor vectors are used whereas in Fig. only 
vectors whose phase trajectories pass through the phase 
tube with length r = 110 are used. It is clearly seen 
that in the first case the vectors z^ and y^ differ from 
each other sufficiently testifying the failure of the pres- 
ence of the functional relation between the interacted sys- 
tem states. But, conversely, for the phase tube with the 
length r = 110 (Fig. ^b) the calculated vectors Zi{t) 
are in the excellent agreement with the vectors yi{t) of 
the second generator that confirms the results obtained 
above for unidirectionally coupled Rossler systems. So, 
in the systems with a mutual type of coupling the vector 
states of the interacting systems are related with each 
other by the functional. 

In conclusion, we have reported that the concept of 
generalized synchronization (except for the LS and CS 
regimes) needs refining, since the state vectors of the 
interacting chaotic systems are related with each other 
by the functional^ but not the functional relation as it 
was assumed until now. Although in the Report the sys- 
tems with a small number of degrees of freedom have 
been considered, the developed formalism can be also 
extended to the systems with the infinite-dimensional 
phase space [25|. Fortunately, this modification of the 
generalized synchronization concept does not discard the 
majority of the obtained hitherto results concerning GS. 
At the same time, this refinement has a fundamental sig- 
nificance from the point of view of the understanding of 
the core mechanisms of the considered phenomena and is 
supposed to give a strong potential for new approaches 
and applications dealing with the nonlinear systems. 
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